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Abstract 



We consider the standard model with two Higgs doublets with the most 
general Yukawa coupling terms ("type III"). In the model, the neutral-Higgs- 
mediated flavor-changing neutral currents (FCNC's) are allowed, but must be 
reasonably suppressed at low energies of probes. It has been known that the 
existing hierarchies of quark masses render this suppression at low energies 
rather natural. On the other hand, the model has been regarded by many as 
unnatural because of the absence of any symmetry that would ensure persis- 
tence of this suppression as the energy of probes increases. The opinion has 
been based on the expectation that the mentioned FCNC's would increase by 
large factors at increasing energies. We performed a numerical analysis of the 
flow of these FCN coupling parameters as governed by the one-loop renormal- 
ization group equations (RGE's), in a simplified case when Yukawa couplings 
of the first quark generation are neglected. The analysis shows a remarkable 
persistence of the mentioned FCNC suppression and thus indicates that the 
model is not unnatural in the RGE sense. Further, we point out two mistakes 
in the Yukawa RGE's of Machacek and Vaughn at one-loop level. 
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I. INTRODUCTION 



The (standard) model with two Higgs doublets whose Yukawa couplings in the quark sec- 
tor have the most general form was apparently first introduced already in 1973 by T. D. Lee 
His main motivation for introducing the model, later also known as the general two- 
Higgs-doublet model (G2HDM) or "type III" 2HDM, was to study new possible CP-violating 
phenomena. Others continued investigating phenomenology of the model along these 
lines. 

Subsequently, Glashow and Weinberg 0] in 1977 stressed that only those models with 
two Higgs doublets whose Yukawa coupling sector possesses specific discrete [or equivalently 
U(l)-type] family symmetries lead to automatic and full suppression of the effective flavor- 
changing neutral (FCN) Yukawa parameters, and ensure this suppression at any energy of 
probes. They pointed out that there are basically two types of such 2HDM's - so called 
"type I" and "type H" models, in which either one Higgs doublet alone is responsible for all 
the quark masses [2HDM(I)], or one Higgs doublet is responsible for all the up-type quark 
masses and the other for the down-type quark masses [2HDM(H)]. This point of Glashow and 
Weinberg apparently had a great impact on the physics community, especially because most 
of the mentioned flavor- changing neutral currents (FCNC's) mediated by neutral HiggsQ 
must be strongly suppressed at low energies of probes (~ Eew) due to the firm experimen- 
tal evidence of FCNC suppression. Consequently, the general 2HDM, which has no such 
automatic and full suppression of the FCN Yukawa coupling parameters, apparently wasn't 
investigated by physicists until the late eighties. 

Since the late eighties, there has been a moderate resurgence of investigation of the 
G2HDM [g- 0. These works investigate low energy phenomena {Er^Eew) as predicted by 
G2HDM's, with most of the FCN Yukawa coupling parameters (at low energies of probes) 
being generally nonzero but reasonably suppressed.^ Refs. [|] investigate predictions of 
the model mainly for CP-violating, and Refs. |]^ mainly for FCNC-violating phenomena. 
The resulting amplitudes then include FCN Yukawa coupling parameters at low energies of 
probes. 

The authors Cheng, Sher and Yuan (CSY) |^ offered arguments which render the 
G2HDM reasonably natural from the aspect of low energy physics, thus countering one 
part of the reservations based on the arguments of Glashow and Weinberg CSY basi- 
cally proposed specific ansatze for the Yukawa parameters in the G2HDM at low energies 
of probes, specifically the FCN Yukawa coupling parameters, motivated laregly by the ex- 
isting mass hierarchy of quarks. Therefore, their ansatze are reasonably natural, or, more 
conservatively, not "unnatural." Motivation of their ansatze didn't explicitly involve any 
family symmetries. Moreover, they showed that their ansatze allow the masses of neutral 



^ A more precise expression would be "neutral flavor-changing scalar (Yukawa) coupling," since 
these couplings have no four- vector current structure involving 7^. 

^ Low energy experiments show that those flavor-changing neutral coupling parameters which 
don't involve t quark are suppressed in nature at low energies E^ruq, while for those involving t 
quark there is no experimental evidence yet available. 
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scalars to be as low as ~10^ GeV while still not violating the available (low energy) data on 
suppressed FCNC phenomena. Later on, Antaramian, Hall and Rasin (AHR) P] proposed 
somewhat similar (but not identical) ansatze, which they motivated by their requirement 
that the Yukawa interactions have certain approximate flavor symmetries. The CSY and 
similar ansatze were mainly used by other authors Q and |^ in their investigations of low 
energy phenomenology of the G2HDM. 

We wish to reemphasize that the mentioned ansatze countered only one part of the 
arguments (based on 0) against the G2HDM. The symmetry arguments of Glashow and 
Weinberg |^ didn't just suggest that a natural 2HDM should have a well motivated sup- 
pression of the flavor- changing neutral (FCN) Yukawa parameters at low energies, but that 
these FCN parameters should remain suppressed also when the energy of probes increases. 
Based on this latter point of , a large part of physics community has continued regarding 
the G2HDM as unnatural. The main point against the G2HDM has consisted of the fear, 
or conjecture, that the FCN Yukawa coupling parameters in the G2HDM, even though sup- 
pressed at low energies by reasonably motivated arguments, would behave unnaturally as 
the energy of probes increases. Stated otherwise, it has been expected that at least some 
FCN Yukawa coupling parameters would increase by a large factor or even by orders of 
magnitude at increased energies (well below the Landau pole), due to the absence of ex- 
plicit discrete [or U(l)-type] family symmetries in the Yukawa sector. The absence of such 
symmetries, so the argument, would in general result in a strong "pull-up" effect on the 
small flavor-changing by the much larger flavor- conserving Yukawa coupling parameters as 
the evolving energy increases. In such a case, the model would then generally contain a 
(thus unnatural) fine-tuning: large "bare" FCN Yukawa coupling parameters at high ener- 
gies would have to be fine-tuned in order to obtain at low energies their phenomenologically 
acceptable suppression. 

Therefore, we investigate this question in the present paper, by performing a numeri- 
cal analysis of the one-loop renormalization group equations (RGB's) of the G2HDM. In 
Sec. 2 we present the model and write down conditions for the suppression of FCN Yukawa 
coupling parameters at low energies (CSY ansatz). In Sec. 3 we write down the one- loop 
RGB's for the Yukawa coupling parameters in the G2HDM in a specific form convenient for 
numerical analyses. A short derivation of the RGB's is given in the Appendix. Sec. 3 con- 
tains also comparison with the existing literature on RGB's. In Sec. 4 we then numerically 
investigate the RGB evolution of the Yukawa coupling parameters for quarks, neglecting for 
simplicity the Yukawa parameters of the light first generation of quarks. We present the 
resulting evolutions of the FCN Yukawa coupling parameters for various low energy ansatze, 
i.e., essentially for variations of the CSY ansatze. We also observe some other interesting 
properties of the presented evolution. Sec. 5 contains a summary and conclusions. 

II. THE MODEL AND LOW ENERGY ANSATZE 

Yukawa interactions for quarks in the G2HDM in any SU(2)L-basis have the most general 
form 
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+U^^\lf¥^^)u%^ + u\fQf¥^^)u%^ + h.c.} + {m-teiYiis} . (1) 

The tildes above the Yukawa couphng parameters and above the quark fields mean that 
these quantities are in an arbitrary SU(2)L-basis (i.e., weak basis, not the mass basis). 
The superscript {E) at the Lagrangian density means that the theory has a finite effective 
energy cutoff E, and the reference to this evolution energy E was omitted at the fields and 
at the Yukawa coupling parameters in order to have simpler notation {E ~ 10^ GeV for 
renormalized quantities). The following notations are used: 



! .„ . (2) 





In (|^), V [= v{Ey\ is the usual VEV needed for the electroweak symmetry breaking, i.e., 
v{Ee^) «i246 GeV. The phase difference ri = i]2 — r]i between the two VEV's in (^) may be 
nonzero; it represents CP violation originating at low energies from the scalar 2HD sector 
(cf. P). The leptonic sector will be ignored throughout. 

We note that the popular "type I" and "type 11" models are special cases (subsets) 
of this framework, with some of the Yukawa matrices being exactly zero: 
[2HDM(I)]; f/(i) = D(2)^o [2HDM(II)]. In these two special models, suggested by Glashow 
and Weinberg Q], the fiavor-changing neutral (FCN) Yukawa coupling parameters are exactly 
zero. This is so because one of the two nonzero Yukawa matrices is proportional to the mass 
matrix of the up-type quarks, and the other to the mass matrix of the down-type quarks. 
Since fiavors refer to the physical quarks, and the quark mass matrices in the physical 
(mass) basis are diagonal by definition, the FCN Yukawa coupling parameters (i.e., off- 
diagonal elements) are zero. Moreover, this is true at any energy of probes (cutoff energy) 
E'. Stated otherwise, when the original cutoff E is changed to E', no loop-induced (ln£" 
cutoff-dependent) FCN Yukawa coupling parameters appear, i.e., the original form of the 
Lagrangian is preserved under the change of the cutoff. This can be formulated also in 
terms of explicit U(l)-type family symmetries governing the Yukawa part of the Lagrangian 
density in the 2HDM(I) and 2HDM(II). These symmetries ensure that, in the course of the 
change of the cutoff (i.e., evolution energy or the energy of probes), the original form of 
the Yukawa Lagrangian density is preserved. In the 2HDM(II) the symmetry tranformation 
is: d^^^ — >-e'"(i^'', — >• e~'"$*^^-' (j = 1,2,3), the other fields remaining unchanged; in the 
2HDM(I): <f^^e"^d^^, u^^ ^q-'^'u^^ , $(2) _,e-i"<|.(2). 

In contrast to the 2HDM(I) and 2HDM(II), the G2HDM has no explicit family symmetry 
enforcing the complete suppression of the FCN Yukawa coupling parameters. There are at 
least two consequences of this fact: 
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1. The FCN Yukawa parameters in the G2HDM are in general nonzero. At low energies 
of probes (£'~i?ew); those FCN Yukawa parameters which don't involve the top quark 
must be given quite small values (not necessarily zero) for phenomenological viability 
of the model. 

2. Even if the FCN Yukawa parameters are zero at some low energy of probes, they 
become in general nonzero at higher energies. If some FCN Yukawa parameters are 
small (but nonzero) at small energies, there exists in principle the possibility that they 
increase by a large factor, or even by orders of magnitude, when the energy of probes 
increases (but remains a safe distance away from the Landau pole) - due to the absence 
of an explicit protective family symmetry. 

The Lagrangian density (0)0 can be written in a form more convenient for consideration of 
the FCN Yukawa coupling parameters, by redefining the scalar isodoublets in the following 
way: 

$'(1) = (cos/5)$(^) + (sin^)e-'''$(^) , 
$'(2) = -(sin/3)$(i) + (cos/3)e-^^<l>(2) , 

where: rj = rj2 — rji ; tan /5 = — ^ cos /3 = — , 

f 1 V 

Therefore, the VEV's of the redefined scalar isodoublets are 

e-.,.(4'...>.^-L(«). ^ijQ. (8) 

The isodoublet is therefore responsible for the masses of the quarks, and couplings of 
to the quarks lead to the FCN Yukawa couplings, as will be seen below. The original 
Yukawa Lagrangian density (|I]) of the G2HDM can then be rewritten in terms of these 
redefined scalar fields as 

41U = - E {GSf (gf $'(^))Ji^-) + Gf (gf I.'«)^i^) +h.c.} 

- E {a,(# $'(^))rfl^ + t/.,(gf $'(^))«g^ + h.c.} , (9) 

where the Yukawa matrices G*^^) and G^^^ are rescaled mass matrices, and U and D the cor- 
responding "complementary" Yukawa matrices, in an (arbitrary) SU(2)L-basis (weak basis) 

G'(^) = ^M^^^jv = (cos/?)f/(^) + (sin/?)e-^''f/(2) , 

G'(^) = V2M^^^ /v = (cos f3)D^^^ + (sin /3)e+"'D(2) ; (10) 
U = -(sin/5)t/(^) + (cos/5)e-'''t/(2) ^ 

D = -{smp)D^'^ + (cos/?)e+'''L'(2) . (11) 



sin P 



V2 
V 



(6) 
(7) 



^ Throughout this Section we omit, for simpler notation, reference to the evolution (cutoff) energy 
E at the quark fields, at the scalar fields and their VEV's, and at the Yukawa coupling parameters. 
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By a biunitary transformation involving unitary matrices , , and V^, the Yukawa 
parameters can be expressed in the mass basis of the quarks, where the (rescaled) mass 
matrices G'^^^ and G^^^ are diagonal and real 



(u) 



V2 



id) 



U 



(12) 
(13) 



VFul , ur = V^ur , d 



VFdL 



d 



R 



V^d 



R 



(14) 



The absence of tildes above the Yukawa coupling parameters and above the quark fields 
means that these quantities are in the quark mass basis (at a given evolution energy E). 
Lagrangian density (|^) can be written now in the quark mass basis. The "neutral current" 
part of the Lagrangian density in the quark mass basis is 



iE)n.c. 
G2HDM 



+Gr^5^4^(0f-i0f)+h.c.} 



V ^ ij=l 



/.'{2)^ 



+ h.c.} . 



(15) 



On the other hand, the "charged current" part of the Lagrangian density in the quark mass 
basis is 



C 



{E)c.c. 
G2HDM 



-(ytG(^)),^.df«0)(^;(i)_i^Ji))+h.c.} 
-4^ E {(VDhufd 



)- 



V2 

-(Vtf/),^J««(^)(0f)_i0f))+h.c.}. 
Here, we denoted by V the Cabibbo-Kobayashi-Maskawa (CKM) matrix 



V = V( 



CKM 



(16) 



(17) 



We see from ([T5|) that the U and D matrices, as defined by (|TT]) and (|T2D-(|T^) through the 
original Yukawa matrices f/^-^-* and D^^^ of the G2HDM Lagrangian density ([1|), allow the 
model to possess in general scalar-mediated FCNC's. Namely, in the quark mass basis only 
the (rescaled) quark mass matrices G^^) and G^^^ of (|12|)-(|13|) [cf. also (|10|)] are diagonal, 
but the matrices U and D in this general framework are in general not diagonal. The 
off-diagonal elements of the matrices U and D are the FCN Yukawa coupling parameters 
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t ) + iu%uM\<pr+^<pr)} ■ (18) 

It should be noted that the original four Yukawa matrices U^^^ and D^^^ (j = 1, 2) in an 
SU(2)L-basis are already somewhat constrained by the requirement that (at low energy) 
the squares o^ M^^^ and M^^^ are diagonalized by unitary transformations involving such 
unitary matrices Vj^ and V^, respectively, which are related to each other by V^V^^ = V. 
Here, V is the CKM matrix which is, for any specific chosen phase convention, more or less 
known at low energies. 

In order to have at low evolution energies (i?~i?cw) a phenomenologically viable suppres- 
sion of the scalar-mediated FCNC's, the authors Cheng, Sher and Yuan (CSY) |^ basically 
argued that the elements of the U and D matrices (in the quark mass basis and at low 
evolution energies E) should have the form: 

U,{E) = ^if^f¥^' , D,iE) = ^^^rr^mf , (19) 

where dj^df-l for E ~ . (20) 

This form is in general phenomenologically acceptable. It is strongly motivated by the actual 
mass hierarchies of the quarks. At least for the diagonal elements, it is suggested by the 
requirement that (at a given low energy ~ -Egw) there be no fine-tuning on the right of 
Eqs. (p!0|)-(p!l|) when these equations are written in the quark mass basis (i.e., no tildes over 
the matrices). For definiteness, consider the up- type sector. The diagonal elements Ujj' are 
in general ^mf^ /v unless fine-tuning is involved on the right of (0). Consequently, also 
Ujj ^ m^j^^ /v unless fine-tuning is involved on the right of (|r^). This consideration further 
suggests (but not necessarily implies) that the off-diagonal elements Uj]^ and Ujk have values 

between those of the corresponding diagonal elements Uj^j ~m^"^/f and ~''^fc'V'^5 
example roughly the geometrical mean of those, leading thus to the CSY ansatz ([T9|)-(pOD.p| 
Therefore, this (CSY) form is considered to be reasonably natural. 

From the CSY ansatz (p!9D-(pOD we see that the FCN Yukawa vertices involving the heavy 
top quark are the only ones that are not strongly suppressed (at low evolution energies). 
As mentioned in the Introduction, scalar-exchange-mediated FCNC processes involving the 
top quark vertices (not loops with top quarks) are not constrained by present experiments. 
Later in Section 4 we will use low energy conditions (|T^-(^Up for a numerical investigation 
of the RGE flow of the FCN Yukawa coupling parameters. 



Strictly speaking, the following "squares": 

^ For the complete suppression of FCN Yukawa couplings Ujk = (for j ^ k) we would then need 
fine-tuning on the right of ([Tl|). 
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III. ONE-LOOP RGB'S IN A CONVENIENT PARAMETRIZATION 



In the Appendix we outlined a derivation of the relevant set of one-loop RGB's for the 
scalar fields and their VEV's (lAloD - dAl^) , for the quark fields ( [AlBD , (|A^ -( [A:2ID , and 
for the Yukawa matrices f/^^^ and D^'^) (^]2|)-(|X]2|). One of the reasons for performing an 
independent derivation is that we consider the method of finite cutoffs , which was used in 
the derivation, as physically very intuitive. This contrasts with other methods often applied 
in the literature, which are, however, usually mathematically more efficient at two-loop and 
higher-loop levels. Another reason is that there is a certain disagreement between the results 
on one- loop beta functions derived for a general (semi) simple gauge group G in various parts 
of literature - see comparisons and the discussion toward the end of this Section. 

We can rewrite all the RGE's derived in the Appendix, now in a a more convenient set of 
parameters. These are: the VEV parameters v = yjvf + V2, tan/5 = V2/V1 and r] = r]2 — rji 

[cf. (^], and matrices G^'^\ U and D [cf. (p!OD, (|1T])] - this representation is more 

convenient for discerning the running of the FCN Yukawa coupling parameters. Applying 
lengthy, but straightforward, algebra to the hitherto obtained RGE's then results in the 
following RGE's in terms of the mentioned set of parameters: 



167r' 



d\nE 



3 Q 9 o 
2^?i + 



(21) 



IGvr 



levr 



, (i(tan/3) 
d\nE 



Nr. 



2 cos2 (3 



d\nE isin(2/5) 



(22) 
(23) 



167r^ 



d 



d\nE 



(U) = nA 2Tr [f/G^^)^ + G^^^D^] G^^^ + Tr \lJlJ^ + IDD^ 



U 



+^iVe(cot/3)[/Tr + f/G(^)t + ^^^^^t - DG^^)^ 



+ < 



1 



-2bb^U - 2G^^^b^G^^'^ - Auil 



(24) 



167r' 



d 



{b) = nJ 2Tr [L)G(^)^ + G(^)f/tl G^^) + Tr + bb^ 



d\nE 



D 



+< 



1 r.- ~ 



tjiji + Db^ + G(^)G(^)t + aiD)QiD)^ 



D + D 



-2UU^b - 2G^^'^WG^^'^ - Aob 



(25) 
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d 



dlnE 



+ 2 



(26) 



167r^ 
1 



1 

+ 2 



+G(^) f^tD + G^^^tG'C^)! _ 2f/G(^)t^ _ 2G(^)G(^)tG'(^) _ An&''\ 



(27) 



Equation (PT| ) is in the Landau gauge, while the other RGE's (P2|)-(P7|) are gauge indepen- 
dent. 

For a general (semi) simple gauge group G, RGE's for various parameters have been 
derived at one-loop level [|10|- and at two-loop level [Q- [jl^. First we should note that 
these groups of authors are using conventions which, particularly as to the fermionic sector, 
differ from each other. 

Cheng, Eichten and Li were using the usual four-component Dirac spinors for quarks. 
While they allowed an arbitrary number of (real) scalar degrees of freedom, their one- 
loop RGE's for the Yukawa coupling parameters are directly applicable only when these 
parameters are real (noncomplex). 



The one-loop RGE results of Vaughn [|TT| and one- and two-loop results of Machacek 
and Vaughn [^] were written for the general case of complex Yukawa coupling parameters. 
Their scalar fields 0a were real, and for fermions (quarks) they were using two-component 
spinor fields ipj as defined by Sikivie and Giirsey |T^. The Yukawa Lagrangian density was 
written in [ITT 



121 in the form 



Cy. = -Y%ipJia^^k<Pa + h.c. 



(28) 



where (pa are the real scalar degrees of freedom, and for fermions (quarks) the two-component 
left-handed spinor fields ipj are used as defined in Ref. [T^; is the second Pauli matrix. 
The double sum over spinor degrees of freedom in (^) is running over j,k = 1, . . . , 2n, 
where n is the number of fermion (quark) flavors. The usual four-component Dirac spinor 
fields '^^^^ in the chiral basis (i.e, the basis of Ref. |]TB|]) are then, according to Sikivie and 
Giirsey 



^ Sikivie and Giirsey ||T5| use the notation Fa, Fb for ipj, ipk, respectively, and (— F* 
Y", (pa, respectively. 



and (hi for 
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^0) 



(J 



,n) . 



(29) 



4jj 

If we carefully rewrite our RGE's (from the Appendix) for the scalar and quark fields, and 
for the Yukawa parameters, in terms of the spinor fields of |T^ and of real scalar field 
components, we notice several differences when compared with one-loop results of |[Tl|]- ||12|| . 
We can deduce from our RGE's for the Yukawa coupling parameters that RG Equation (3.4) 
of ^ (second entry), or equivalently, RG Eqs. (2.2)-(2.3) of H 
in their language should read: 



for Y" Yukawa matrices 



(4vr) 



dlnE 



= (47r)2/5" 



1-1. 



1-1. 



Y^Y^^Y'^+ Y'^Y^^Y^ _l_gY^Y'^^Y''-|- 



+4/tY^Tr (^Y^tY»+Y"tY^) - 3g^ {C2(F), Y"} 



(30) 



where k = 1/2. Stated otherwise, the cubic Yukawa terms on the right of this RGE are 
effectively those given in [|T^, but multiplied^ by factor 4, and the trace there is replaced 
now by its real (symmetric) part. Similar differences arise when comparing our RGE's for 
scalar and quark fields with those of []12[. Instead of Eq. (3.7) of [|12[ (first entry), we get 

4KTr fY"Y^^ + Y'^^Y^^ 



(47r)^7:.|i-i. 



g\2+a)C2iS)5at 



(31) 



where again k = 1/2 and 7^* is defined via: d(f)a/d\n E = —'j^f^(j)h- Here, (pa are the real scalar 
fields, and a in (|3T|) is the gauge parameter (a = l — ,^ = 1 in the Landau gauge). Instead of 
Eq. (4.5) of PI (first entry), we get 



(477)^7 



1-1. 



2Y'^tY« + /C2(F)(l- 



a) 



(32) 



where 7^ is defined via: dipj/dlnE 



spinors as defined in |T5] 
The authors of 



Here, ipi are left-handed two-component 



13 



on the other hand, worked with Majorana fermions, using back- 
ground field method. The Dirac fermions can then be expressed as sums of two Majorana 
fermions. Their one-loop beta functions for the Yukawa coupling parameters of the (real) 
scalars with the (Majorana) fermions, can be re-expressed in the notation with left-handed 
two-component spinors ipj as introduced by |]15| and used by [|lT|- [1^. After a somewhat 
lengthy algebra, it can be shown that the one-loop results of fl^ lead precisely to formula 
(BO). Therefore, we finally conclude that our one- loop RGE formulas for the Yukawa cou- 



pling parameters, derived in the Appendix and rewritten in Eq. (^) in the language of [|15 
are not in agreement with those of Vaughn |11| and of Machacek and Vaughn , and are in 



agreement with the results of Jack and Osborn |T^. Moreover, the latter authors emphasize 
that their RGE results agree with those of van Damme ||14||.p| 



^ If Machacek and Vaughn had introduced in the Lagrangian density an additional factor of (1/2) 
in front of the sum (|2^ ) ( - but they didn't), the factor 4 in the cubic Yukawa terms on the right 
of (^) wouldn't have occurred. 

^ Here we also mention that Fischler and Oliensis have derived RGE's for Yukawa coupling 
parameters of the minimal SM at two-loop level. 
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For several reasons, we considered it instructive to perform an independent derivation 
of the one-loop RGB's for the scalar and quark fields and for the Yukawa matrices in the 
discussed general 2HDM. One reason is that the one-loop results of [|1T| and 
agree entirely with those by other authors 



12 do not 
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Another reason is that the existing 
works on the one- and two-loop RGB's for general (semi)simple gauge groups G use various 
conventions for the fermionic fields, and are usually written in a language difficult for non- 
specialists in the method used. The third reason is that these works do not apply the method 
of finite cutoffs which we consider especially appealing and physically intuitive - although, 
at two-loops, probably not the most efficient one. With our independent cross-check we are 
confident that the one- loop results of [|13| are correct. 



IV. NUMERICAL EXAMPLES OF EVOLUTION 

Here we present a few simple but hopefully typical examples of the RGB evolution of 
parameters in the G2HDM. Some preliminary numerical results were presented by us in 



Ref. [18|. For simplicity, we assume: 



• there is no CP violation - all original four Yukawa matrices U'^^\ D^^^ are real, and 
the VBV phase difference r] is zero; 

• the Yukawa parameters of the first quark generation as well as those of the leptonic 
sector are neglected (the quark Yukawa mass matrices are therefore 2x2). 

For the boundary conditions to the RGB's, at the evolution energy E = Mz, we first take 
the CSY ansatz (|19D-(|0D, with ^^f = 1 = df or ^^f = 2 = i\f for alH, j = 1, 2. We stress 
that i = l refers now to the second quark family {c,s), and i = 2 to the third family {t,h). 
For the (2x2) orthogonal CKM mixing matrix V we take Vi2{Mz) = QMb = -V2i{Mz). The 
values of other parameters at E = Mz are chosen to be: 
tan/3 = 1.0; t; = y^t;2+t,2 = 246.22 GeV; 
03 = 0.118, 02 = 0.332, ai = 0.101; 

me = 0.77 GeV, m, = 0.11 GeV, m6 = 3.2 GeV, and mt = 171.5 GeV. 

The latter quark mass values correspond to: mdrnc) ~ 1.3 GeV, ms(lGeV) ~ 0.2 GeV, 
mb{mb) ^4.3 GeV, and m^^^^' ^ 174 GeV [mt{mt) ~ 166 GeV]. For a^lE) we used two-loop 
evolution formulas, with threshold effect at E^rrii^^^' taken into account; for aj{E) {j = 1,2) 
we used one-loop evolution formulas. 

The described simplified framework resulted in 18 coupled RGB's [for 18 real parameters: 
f^, tan/3, Uij, Dij, G\'^\ G^f^], with the mentioned boundary conditions at E = Mz- The 
system of RGB's was solved numerically, using Runge-Kutta subroutines with adaptive 
stepsize control (given in |l^). The numerical results were cross-checked in several ways, 
including the following: FORTRAN programs for the RGB evolution and for the biunitary 
transformations were constructed independently by two of the authors (S.S.H. and G.C.), 
and they yielded identical numerical results presented in this Section. 

The results for the FCN Yukawa parameter ratios Xij{E) / Xij{Mz) {X = U,D; i^j) are 
given for the case = l = ^lf in Fig. |I]. From the Fi gure we immediately notice that the 
FCN coupling parameters are remarkably stable as the energy of probes increases. Bven 
those FCN Yukawa coupling parameters which involve t quark remain quite stable. Only 
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FIG. 1. J^tlN Yukawa parameter ratios Ufff^ 



D,j{E)/D,j{Mz) {i ^ j) in tHe G2HDM as 



the Euclidean energy of probes E increases. These parameters are in the quark mass basis. At E = Mz, 
CSY ansatz was taken with ^^-"^ —^i? — 1 (for aU ij — 1,2). 



very close to the top-quark-dominated Landau pole (£'poie~0.84 ■ 10^^ GeV)0 the coupling 
parameters start increasing substantially. For example, in the down-type FCN sector (6- 
c) the corresponding ratio D2i{E) / D2i{Mz) acquires its double initial value (i.e., value 2) 
at ~ 0.7i?poie, which is very near the (Landau) pole. About the same holds also for 
U2i{E) / U2i{Mz) . For the ratio Dui^E) / D^i^Mz) the corresponding energy is even closer to 
Epoie- For the t-quark-dominated Uui^E) /Uui^Mz) it is somewhat lower. 

In Fig. 0, evolution of the same FCN ratios is depicted for the case of the low energy 



CSY parameters ^f"^ = 2 = ^),"''. The t-quark-dominated Landau pole is now substantially 



5ij 



5 



lower (i?poie~10 TeV), but the behavior of the FCN ratios remains qualitatively the same. 
Moreover, when some of the CSY parameters are varied, the stability of the FCN ratio 
persists, and the Landau pole is influenced almost entirely by the t-quark-dominated CSY 
parameter ^22'' • We also looked into cases when the CSY ansatz is effectively abandoned. 
If we suppress the up-type off-diagonal element at -E = Mz drastically, for example by 
taking ^[2 =^21 ~ 0.0516 [corresponding to Uij{Mz) = Dij{Mz) for i^j] and all other 
parameters equal to 1, we obtain results depicted in Fig. ^, which are very close to those of 
Fig. |. 

From all these Figures we conclude that the FCN Yukawa coupling parameters in the 
general 2HDM show remarkable stability when the (Euclidean) energy of probes increases. 
This stability persists up to the energy regions which are, on the logarithmic scale, quite 
close to the top-quark-dominated Landau pole. The general 2HDM appears to possess this 
stability even when the off-diagonal low energy parameters \C,ij\ 
smaller than 1 (while the diagonal parameters are ^jj~l). 



m- 



have values much 



^ The value of ii'poie is strongly dependent on the given value of parameter ^, as shown later in 
Fig. |. 
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FIG. 2. Same as in Fig. but for the choice ^'^f^^'^f = 2 (for aU i,j = 1,2) 



We can also compare the RGE evolution of the Yukawa parameters in the G2HDM with 
those of the MSM, 2HDM(II) and 2HDM(I). The authors of |lj proposed that (heavy) 
quark mass and Higgs masses in the MSM could be determined by the infrared fixed points 
of the RGE's. These questions were numerically investigated also in different variants of 
the 2HDM(II) and 2HDM(I) |2|]. The authors of the latter work found out that relatively 
unambiguous predictions can be made only if there is a heavy quark generation and the 
(heavy) quarks couple to both Higgs doublets. It is interesting to note that in the G2HDM, 
the heaviest {t) quark also has an infrared fixed point behavior, as suggested from Figs. |^-|^. 
This is further suggested from Figs. |^-|^ which represent evolution of the ratios X{E) / X{Mz) 
involving the flavor-nonchanging neutral Yukawa parameters {X = Ujj or Djj) or the "mass" 
Yukawa parameters {X = G^^^ or G^^^). By the infrared fixed point behavior we mean that, 
for a given approximate Landau pole energy Epo\e, we have very weak sensitivity of X{Mz) 
on the otherwise large value X(ii^poie) ~ 1. The reason for this similarity with the MSM and 
2HDM(II) and 2HDM(I) lies probably in the conjunction of the facts that the CSY ansatz 
(P^)-(P(]|) implies dominance of the t-related Yukawa coupling parameters and that the QCD 
contribution to the evolution of these parameters has the sign which is in general opposite 
to that of the (t-related) Yukawa parameter contributions. Inspecting the work p2[ further, 
we note that it would also be interesting to investigate the RGE behavior of the quartic 
Higgs coupling parameters in the G2HDM. This would tell us when these parameters have 
infrared fixed point behavior und thus when the physical Higgs boson mass spectrum can 
be determined. 

In this context, we mention that the idea of RGE fixed points had been introduced earlier 



by Chang |23]. He, and subsequently others |24|, investigated connection of RGE fixed points 



with asymptotic freedom in massive gauge theories. Cabibbo et al. were apparently the 
first to investigate mass constraints in the minimal SM by imposing boundary conditions on 
(perturbative) RGE's at the unification energy of grand unified theories [SU(5) or 0(10)]. A 
somewhat related analysis was performed later by the authors of Ref. p6[ who used SU(5) 
fixed-point conditions. 
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FIG. 3. Same as in Fig. but for the choice ^Y^' =0.05163 (other ^^'s are 1) 



It should be stressed that the presented numerical results are independent of the chosen 
value of the VEV ratio, tan/?, at E = Mz- This is connected with our choice of the CSY 
boundary conditions (p!9|)- (|20|) at E = Mz for the Yukawa matrices in the quark mass basis 
(all ^jj's taken real) and the reality of the chosen CKM matrix at E = Mz- These boundary 
conditions result in real and /3-independent Yukawa matrices U, D, G^^^ in a weak 

[SU(2)l] basisQ at E = Mz- The RGB's (g)-(|23) then imply that these matrices remain 
real and independent of (3 at any evolution energy E, and that also their counterparts U, 
D, G(^) and G^^) in the quark mass basis, as well as the CKM matrix V, remain real and 
independent of (3 at any energy E. Stated otherwise, if there is /3-independence and no CP 
violation (neither in the original Yukawa matrices nor in the scalar sector) at a low energy 
{E = Mz), then these properties persist at all higher energies of evolution.[^ 

This feature is in stark contrast with the situation in the 2HDM(II) where the Yukawa 
matrices strongly depend on (3 already at low energies - e.g., gt{Mz) = mt{Mz)V2/vu = 
mt{Mz)V2/[v sin (3 {Mz)]. Also the location of the Landau pole in the 2HDM(II) then 
crucially depends on j3{Mz) - smaller (3{Mz) implies larger gt{Mz) and hence a drastically 
lower Landau pole. 



10 We chose at E = Mz the following weak basis: U = U, G^^^ =G^^\ D = VD, = VG^^\ 
where V is the CKM matrix (at E = Mz)- According to relations (|lO|)-([Tl|), the reaUty of the 
Yukawa matrices U , D, and G^^^ at the low energy E = Mz would follow, for example, from: 
the requirement of no CP violation in the Yukawa sector (i.e., the original Yukawa matrices U^^^ 
and D^^^ are all real) together with the requirement of no CP violation in the scalar sector (i.e., 
the VEV phase difference r/ = 0) at that low energy. 

11 CP conservation in the pure scalar sector at a low energy E = Mz (i.e., f? = 0) also persists 
then at all higher energies of evolution, since drj/d\nE = t) by the reality of the Yukawa matrices, 
according to RGE (H). 
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FIG. 4. Same as in Fig. but for the neutral current Yukawa coupling parameters Ujj and Djj [j = 1,2) 
which don't change flavor. 



On the other hand, the G2HDM treats the up-type and the down-type sectors of quarks 
(the two VEV's Vi and V2) non-discriminatorily. Therefore, it should be expected that 
any reasonable boundary conditions for Yukawa coupling parameters at low energies should 
also be independent of (3 in such frameworks, and this independence then persists to a 
large degree also at higher energies. Also the locations of the Landau poles (i.e., of the 
approximate scales of the onset of new physics) should then be expected to be largely j3- 
independent. In this sense, the G2HDM has more similarity to the minimal SM (MSM) 
than to the 2HDM(II). The persistence of complete /3-independence of the Yukawa coupling 
parameters at high energies and of the Landau poles, however, can then be "perturbed" 
by CP violation - because RGE's (p^)- (|27|) are somewhat /3-dependent when the Yukawa 
matrices U, etc., are not real. 

In addition to the connection between (low energy) CP violation and /^-dependence of 
high energy results, there is yet another feature that distinguishes the G2HDM from the 
MSM - the Landau pole of the G2HDM is in general much lower than that of the MSM. 
We can see that in the following way: let us consider that only the Yukawa parameters 
connected with the top quark degree of freedom are substantial, i.e., G22'' = 9t ^ ^ and 
U22 = g't^l. We have: gt{E) =mt{E)y/2/v{E), as in the MSM, and g'tiE) is an additional 
large Yukawa parameter - both crucially influence location of the Landau pole. Inspecting 
RGE's ( P^ and ( PB| ) for this special approximation of two variables gt and g'^, we see that 
RGE for gt is similar to that in the MSM, but with an additional large positive term on the 
right: {3/2){g^)'^gt. RGE for g[ has a similar structure as RGE for gt, but with substantially 
larger coefficients at the positive terms on the right. As a result, g'tiE) is in general larger 
than gt{E). Our specific numerical example ^fj'' = ^if = 1 shows that g'tiE) is on average 
(average over the whole evolution energy range) almost twice as large as gt{E). If we then 
simply replace in the mentioned additional term {3/2){g't)'^gt the parameter (g't)"^ by 3.5^^^, 
we obtain from the resulting "modified" MSM RGE for gt a value for the Landau pole in 
the region of 10^^-10^^ GeV, which is roughly in agreement with the actual value of the 
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logio(Energy/GeV) 

FIG. 5. Same as in Fig. |l|, but for the "mass" Yukawa parameters G^'^^' 



and G 



[D] 



^..^ {j = 1,2) instead. 

Since G'^^^E) and G'^^^E) matrices are diagonal by definition (quark mass basis), these neutral current 
Yukawa matrices have zero FCN components automatically. 



Landau pole of our numerical example -E'poie~0.84 ■ 10^^ GeV. And this value is much lower 
than -Epoie in the MSM which is above the Planck scale. As already mentioned, the value 
of -Epoie is largely influenced by the value of the top-quark-dominated parameter ^22'*- Of 



course, when we allow the (^22'') parameters of the CSY ansatz (|19 D -(^0 D at E = Mz to 



deviate from 1, we obtain larger log(i?poie) for smaller ^jj, and smaller log(£'poie) for larger 
^ij. In Fig. ^we depicted this variation of the Landau pole energy when the CSY low energy 
parameters are varied. 

It is also interesting to note that RGE (|23|) for the evolution of the difference of the 
VEV phases {rf) implies in the G2HDM that 77 can change when the energy of probes 
changes. This would generally occur when we have CP violation in the Yukawa sector (i.e., 
complex Yukawa matrices). Even when 7^ = at some low energy, it may become nonzero at 
some higher energy due to the CP violation in the Yukawa sector. This contrasts with the 
2IIDM(II) or 2HDM(I) where the right side of (^Sf) is zero always and thus the CP violation 
in the Yukawa sector doesn't influence 77. 



V. SUMMARY AND CONCLUSIONS 

We performed numerical analysis of the one-loop RGE's in the general two-Higgs-doublet 
model (G2IIDM). In the analysis, we neglected the Yukawa coupling parameters of the 
light first quark generation, as well as the contributions of the leptonic sector. At low 
energies of probes, we first adopted the CSY ansatz (|19|)-(^0|) which is largely motivated 
by the existing quark mass hierarchies. We found out that the flavor-changing neutral 
(FCN) Yukawa parameters remain remarkably stable when the energy of probes increases 
all the way to the vicinity of the (t-quark-dominated) Landau pole. This conclusion survives 
even when the CSY ansatz is effectively abandoned, i.e., when the off-diagonal low energy 
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1 1.5 2 2.5 

FIG. 6. Variation of the Landau pole energy when the low energy parameters ^l"'' —^^j ^=£^ of the CSY 
ansatz (pj|)-(pO|) are varied. For ^ = 2.5, the onset scale of new physics is already quite low: Epoic~2 TeV. 

parameters are additionally suppressed: <^ 1 {i^j). This behavior indicates that the 
G2HDM doesn't behave unnaturally with respect to the RGE evolution of the vertices of 
the Higgs-exchanged flavor-changing neutral currents. Since the G2HDM, in contrast to 
the 2HDM(II) and 2HDM(I), has no explicit and exact discrete [or U(l)] family symmetries 
which would ensure persistence of the FCN Yukawa suppression at increasing energies of 
probes, the behavior of FCN Yukawa parameters found numerically in the present paper may 
be somewhat surprising. The general suspicion about the G2HDM in the past had centered 
on the fact that absence of the mentioned family symmetries in the Lagrangian density 
would in general not keep FCN Yukawa parameters suppressed under the RGE evolution 
and would thus render the model unnatural and fraught with flne-tuning of "bare" FCN 
Yukawa parameters. Stated otherwise, the RGB's of the G2HDM would in general allow 
a "pull-up" effect by the diagonal Yukawa parameters on the much smaller off-diagonal 
(FCN) ones. This would increase the values of the latter by a large factor or even by orders 
of magnitude when the energy of probes increases by one or several orders of magnitude. 
Our numerical analysis shows that this doesn't happen, at least as long as the low energy 
parameters of the CSY ansatz (|19D-(^0D saflsfy \ ~ 1 for iy^j and ^jj ~ 1. Perhaps it is 
the latter condition (for the diagonal ^ parameters) which causes the mentioned persistence 
of the FCN Yukawa suppression. The latter condition, together with the known form of the 
CKM matrix, effectively represents an approximate symmetry in which only the third quark 
generation has substantially nonzero Yukawa parameters (and almost no CKM mixing). 
This can be called also approximate flavor democracy. 

Further, the high energy Yukawa coupling parameters in the model have in general little 
dependence on the VEV ratio tan /3 as long as the CP violation is weak. Moreover, we found 
out that the G2HDM has an interesting behavior of the Landau pole energies. They can 
become quite low (~ 1 TeV) already at not very high C,ij parameters (^22'' < 3); as shown 
in Fig. ^ These energies, signaling the breakdown of the perturbative approach in the 
G2HDM, can be interpreted as possible scales of the onset of a new physics and/or a strong 



17 



coupling regime. 
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FCNC - flavor-changing neutral current; 

G2HDM - general two-Higgs-doublet (Standard) model; 

MSM - minimal Standard Model; 

RGE - renormalization group equation; 

VEV - vacuum expectation value; 

IPI - one-particle-irreducible; 

IPR - onc-particlc-rcducible; 

2HDM - two-Higgs-doublet (Standard) Model. 
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Appendix A. ONE-LOOP RGE'S IN THE GENERAL 2HDM 



We outline here a derivation of the one-loop RGE's for the scalar and quark fields and 
for the Yukawa coupling matrices D^'^^ and JJ^'^^ (fc = l,2) in the general two-Higgs doublet 
model (G2HDM) whose Lagrangian density in the Yukawa sector is represented by ([l|). The 
derivation follows the finite cutoff interpretation of RGE's as presented, for example, in 
Ref. [0. We consider it useful to present the derivation because the approach used here is 
physically very intuitive, and because the existing relevant literature on RGE's in general 
(semi) simple Lie gauge groups is written in a rather cryptic manner and not all the works 
agree completely with each other. In Sec. 3 we also compared the results obtained here with 
those implied by the existing literature. 

A.l. One-loop RGE's for the scalar fields 

To obtain evolution of the scalar fields (j)'f^\E) with "cutoff" energy E, the trun- 
cated (one-loop) two-point Green functions — iEj^'^^(p^; ii^^), represented diagrammatically 
in Fig. 0, have to be calculated first. More specifically, it suffices to calculate only their 




q 

FIG. 7. The diagram leading to the two-point Green function —\Y!f^'^\p'^ ; E"^) . Full lines represent 
quark propagators. 

cutoff-dependent parts oc In E"^ which are responsible for effective^ kinetic-energy-type 
terms ^ d,y(f)[''\E)d'^(j)f\E). In the course of the calculations, all the masses m ~ i?ew of 
the relevant particles in the diagrams are ignored. This would correspond to the picture 
with a finite but large ultraviolet energy cutoff E^E^^. Therefore, calculations need not 
be performed in the mass basis of the relevant particles. These particles are regarded as 
effectively massless in the approximation, the transformations between the original bases of 
the relevant fields and their mass bases are unitary, and therefore the (mass-independent 
parts of the) calculated Green functions are the same in both bases. 

Calculation of the Green functions —iTif'j^\p^; E"^) is then straightforward. The relevant 
massless integrals over internal quark-loop momenta q can be carried out in the Euclidean 



For clearer notation, we denote in this Section the evolving (UV cutoff) energy E at the fields 
not superscript, but rather as an argument. 
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metric [q = (— g-^), p = {—ip^,—p^)], where the upper bound in the loop integral is: 
(f<E'^. After rotating back into Minkowski metric {p^^—p^)^ we obtain: 



1. Green functions whose external legs cjyf' and (jy'f^ have the same scalar indices {i=j): 



xTr 



where j = 1,2,3,4 (no running over j); k,i = 1,2; m is an arbitrary but fixed mass 
(m~Eew). 

(k) (i) 

2. Green functions whose external legs (f))j and (f)j, have complementary scalar indices 
0y) = (12), (21), (34), (43): 



X 



xTr 



t/Wf/Wt _ ifWum _ ^(fc)^Wt + ^W£)Wtl (^) . (A.2) 



3. Green functions — iS,-^'' with other indices are zero. 

All these Green functions can be induced alternatively at the tree level by kinetic 
energy terms. For example, in theory with the UV cutoff E, the kinetic energy 
term du(j)f'\E)d''(j)f\E) induces (at the tree level) the two-point Green function value 
-iJ:i';/\p^;E^) = ip2 if 0f)(E) ^ (pfiE), and the value 2ip2 if 0f)(E) = (pfiE). Now, 
following the finite-cutoff interpretation of RGB's as described, for example, in Ref. [^, we 
compare the kinetic energy terms in the theory with the UV cutoff E and in the equivalent 
theory with the slightly different cutoff (E + dE). The two-point Green functions in these 
two equivalent theories must be identical. Imposition of this requirement in the tree plus 
one-loop approximation then leads to the following relation: 

^ E E d4'\E)d^cl^f\E) = ^ E E d^<j>f\E+dE)d'^<j>f\E + dE) + 

^ j=l k=l ^ j=l k=l 

+-^(dlnE^)|E E A,,{E)d.^f\E)d'^<pf{E) + 
•^^^ I j=i k,e=i 

+ E E {-iy'^'BUE)d,<pf\E)d^^fiE)\ , (A.3) 

where: summation in the last sum runs over already mentioned complementary indices 
(jj') = (12), (21), (34), (43); d{\n E^) = \n{E+dEf-\n E^ = 2dE / E; and elements of the real 
symmetric matrices A{E) and B{E) are related to the Green function expressions ( |A. 1|) and 
( |A.2| ), respectively: 
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Ak,{E) = ^Tr [[/«[/Wt + f/Wf/(fc)t + ^W^Wt + ^W^Wt] (^) ^ 
Bm{E) = (-l)^^Tr [t/W?7Wt_f/Wt/Wt_^W^(m + ^W^(fc)tl (^) . 



(A.4) 
(A.5) 



Equation ( |A.^j| ) is described in the following way: the double sum on the left and the first 
double sum on the right represent the kinetic energy terms of the scalars in the formulation 
with UV cutoff E and {E+dE), respectively. The one-loop contributions of Fig. ^ with the 
loop momentum \q\ in the Euclidean energy interval E <\q\<A are already contained in the 
kinetic energy terms of the left effectively at the tree level (A is a large cutoff where the theory 
is presumed to break down). On the other hand, the kinetic energy terms of the (E + dE) 
cutoff formulation [the first double sum on the right of (|A.3|) ] effectively contain, at the tree 
level, the one-loop effects of Fig. |^for the slightly smaller energy interval: (E+dE) < \q\ < A. 
Therefore, the Green function contributions]^ —idTi\j^\p'^;E'^) of Fig. from the loop 



momentum interval E<\q\< {E + dE) had to be included on the right of ([A.3|) - these are 
the terms in the last two double sums there. This is illustrated in Fig. 0. 



p+q 



ij,k,l 



E<lql<E+dE 



A. ' 

- 

P 



a) 




(I) 



= 



FIG. 8. Diagrammatic illustration of the RGE relation (|A.3| ) leading to the evolution of the scalar 



fields. 0f'^ stands for (l>'i'''{E), and dcj) stands for 4){E + dE) ~ (f>{E) is a generic notation for i/i^- ■^'s). The 
cross represents the contribution of the change of the kinetic energy terms originating from the changes dcj) 
of scalar fields. 



In order to find RGE's for the scalar fields 
solution of Eq. ( |A.3| ): 



(E), we make the following ansatz for 



(pj{E + dE) = (Pj{E) + da^^\E)$j{E) + dj3^^\E)$j,{E) 



(A.6) 



where 0j is two-component column made up of cpj and 0^- ; da'^^^E) and d/S'^^^E) are 
infinitesimally small 2x2 matrices, and dp^^^ has zero diagonal elements; and scalar indices 
{jj') are again complementary. Inserting ansatz ([A.6| ) into RGE relation ( A. 3 ), we obtain 
relations 



dag(^) + da%\E) = -N, [d\nE^) Am{E) / (IQir^) , 



(A.7) 
(A.8) 



In principle, these relations alone do not define the elements da'i}{E) and df3'^l){E). However, 
RGE evolution of the isodoublet fields ^^^^E) and should be invariant under the 



More precisely: the corresponding effective kinetic energy terms. 
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exchange of Higgs generation indices l<-i>2, because these two Higgs doublets appear in the 
original Lagrangian density (|l|) in a completely 1^2 symmetric manner. We will see in 
retrospect that this exhange symmetry is respected once we impose the conditions 



da^^i {E) = day> (E) , d(3'^> {E) = d(3l > {E) 



Si) 



(A.9) 



Relations (|A.7|) - (|A.9|) lead to specific expressions for the evolution coefficients da'l}{E) and 
d(3^^}{E). When inserting these coefficients back into ansatz ( |A.6| ), we obtain the one- loop 
RGE's for the evolution of the scalar fields: 



Nr d\nE^^ 



-Tr 



i/(i)f/(i)t + ^(i)^(i)t 



1 

— Tr 
2 



■^(l)fj(2)t _^ ^(2)^(l)t _,_ /)(2)t _,_ £)(2) ^(l)t 



A, 



d2) 



+i(-iy ^Tr [t/«t/(2)t _ ^(2)t;(i)t _ ^(i)^(2)t + ^(2)^(i)t 



^(2) 



(A.IO) 



IGtt^ d (2) 



-Tr 



' c 

It. 

2 



+i(_l)i+iiTr [f/«?7(2)t - m^)uW] _ ^(i)^(2)t + ^(2)^(i)t 



(A.ll) 



where again j' is the scalar index complementary to index j: {jj') = (12), (21), (34), (43). 
These RGE's lead to RGE's for scalar isodoublets 



167r2 



Nr dlnE 



(A.12) 



1=1 



as 



We really see that this set of one-loop RGE's is invariant under the exchange 1 ^ 
required by the form of the Yukawa Lagrangian density (|I]). 

In addition to quark loops, there are also loops of the electroweak gauge bosons con- 
tributing to one-loop two-point Green functions of the scalars. However, since these gauge 
bosons couple to the Higgs isodoublets identically as in the minimal Standard Model (MSM), 
their contributions]^ to the right of RGE's ( |A.1CI| )-( [A.11| ) and ( |A. 12| ) are the same as in the 
MSM.|3 Hence, the full one-loop RGE's for the the scalar isodoublets in the G2HDM, in the 
Landau gauge, are 
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They are gauge dependent. 



For these contributions of EW gauge bosons in the MSM, see for example Arason et al. [27|, 
App. A. However, note that they use for the U(1)y gauge coupHng gi a different, GUT-motivated, 

convention: (5i)Arasonetal. = {5/3){gl)here- 
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d 



d\nE ^ ' ^ 



+ 



(A.13) 



and completely analogous gauge boson contributions have to be added also on the right of 
( |A.1CI| )-( ]A.11|) . These RGB's are simultaneously also RGB's for the corresponding VBV's 

(I) 



d 



dlnE 



Vi 



+ 



(A.14) 



In this paper we don't discuss the question of quadratic cutoff terms A^ which appear in 
the radiative corrections to VBV's in any SM framework. In the MSM, their consideration - 
under the assumption of the top quark dominance of the radiative corrections in the scalar 
sector - leads to severe upper bounds on the ultraviolet cutoff A for a substantial subset of 
values of the bare doublet mass and of the bare scalar self- interaction parameters M^(A) 
and A (A) - cf. Ref. 
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In order to derive one-loop RGB's for the Yukawa matrices U^^^ and D'^^\ the results 
( |A.1CI| )-( [A.13| ) are needed. In addition, RGB's for evolution of the quark fields u^^^ and d^^ji 
are also needed. 



A. 2. One-loop RGE's for the quark fields 

These RGB's can be derived in close analogy with the derivation of the evolution of scalar 
fields of the previous Subappendix. Now, the diagrams (Green functions) of Figs. |^ and ^ 
are replaced by those of Figs. ^ and 10, and the scalar field kinetic energy terms in ( |A.3|) 
are replaced by those of the quark fields. The Green function of Fig. 0, with the incoming 



u 



u 



p q-p p 

FIG. 9. The diagram leading to the two-point Green function — iS(p; i?; u''', w'^^). Dashed and full 
lines represent scalar and quark propagators, respectively. 



U 



and outgoing u^^^ of momentum p, in the framework with UV cutoff E, is 



-iT.{p-E- 



6477' 



In (^]A2{l+j,)j:\mm'^ 



1=1 



+ 



(A.15) 
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The Green function with the incoming c?*^*) and outgoing d^-^^ of momentum p is obtained 
from the above expression by simply exchanging U^^'> ^ D^^^ and IJ^^^"^ ^ D^^'^'^ . The quark 
fields evolve according to the ansatz 

= df,{E)[\^'''^q^'\E),^n , (A.16) 

where dq{E) generically stands for q{E + dE) — q(E) [q^''^ = u^''\ d^'^'^], and subscripts L, R 
denote handedness of the quark fields: gz, = (1—75)5/2, 5^? = (1+75)5/2. In complete analogy 
with the previous Subappendix, we obtain from this ansatz and from the RGE relation 
illustrated in Fig. |10|^°| the following relations for the quark field evolution matrices df^. 




E<lql<E+dE 



FIG. 10. Diagrammatic illustration of the RGE relation leading to the evolution of quark fields. This 
relation means that the two-point Green functions with truncated external quark legs, at one-loop level, are 
the same in the theory with E cutoff and in the theory with E+dE cutoff. Conventions are the same as in 
previous figures. 



dUE)\f^* + dUE)\'^' 



dUE)f* + dUE)f 



2\ 2 



327r2 ^' 
2(rflnE2)^j^(,)^^(,)j .^(E). 



327r2 



k=l 



(A.17) 
(A.18) 



The relations for the df^ evolution matrices of the down-type sector are obtained from the 
above by simple exchanges U^^^ ^ D^^'^ and U^^^'^ ^ D^^^"^ . Hermitean solution to these 
relations is 



dfu{E)'lf^ 
dUE)f 



dUE) 



(R) 



{dlnE 
647r2 

2{d\nE 
647r2 

2{d\nE 
647r2 



2\ 2 



k=l 
2\ 2 



k=l 
2\ 2 



(A. 19) 
(A.20) 
(A.21) 



fc=i 



The results ( [A.19| )-( [A.21| ), in conjuction with ( [A.16| ), represent one-loop RGE's for evolution 
of the quark fields in the G2HDM, but without the gauge boson contributions. The latter 
contributions are the same as in the MSM and can be included in ( |A.19| )-( [A.21| ). 
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RGE relation represented by Fig. IC is analogous to relation ( [A.3| ) represented by Fig 
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A. 3. One-loop RGE's for the Yukawa coupling matrices 



To derive these RGE's, we need, in addition to results of the two previous Subappendices, 



also another Green function. It is represented by the diagram of Fig. O. When the external 




FIG. 11. One-particle-irreducible (IPI) diagram contributing to the evolution of the Yukawa coupling 
parameters. Conventions are the same is in previous figures. 



legs there are u^"^' (incoming, with momentum fc), u^^' (outgoing, with momentum p+fc), and 
03^^ [or 04^^], it turns out that only the diagram with the charged scalar exchange contributes, 
and the resulting truncated three-point Green function, in the framework with the UV cutoff 
is 




327r2v^ 




1+75) 



+ (1-75) 



(A.22) 



The corresponding Green function with the down-type quark external legs is obtained from 
the above by the exchanges and U^'^^^ -i^ D^^^"^ . 

Now, the one-loop RGE's for the Yukawa matrices are obtained in analogy with the 
reasoning leading, in the case of two-point scalar Green functions, to the RGE relation ( [A. 3 ) 
[cf. Fig. ||. It is straightforward to check that the contribution of the quark loops in the scalar 
external leg cancel the contributions coming from the renormalizations of the scalar fields in 
the kinetic energy terms of the scalars - this is illustrated in Fig. O. Furthermore, it can be 




FIG. 12. Cancellation of contributions from the quark loop (one-particle-reducible - IPR) with those 
of the scalar field renormalizations in the kinetic energy term of the scalars, for the energy cutoff interval 
{E,E+dE). 

checked that the contributions of the scalar exchanges on the external quark legs cancel the 
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contributions coming from the renormalizations of the quark fields in the kinetic energy terms 
of the quarks - illustrated in Fig. 13, All in all, the IPR one-loop contributions are canceled 







= 



E<lql<E+dE 

FIG. 13. Cancellation of contributions from the scalar exchange on the quark legs (IPR) with those 
of the quark field renormalizations in the kinetic energy term of the quarks, for the energy cutoff interval 
{E,E+dE). 



by the contributions of field renormalizations in the kinetic energy terms. Therefore, the only 
one-loop terms contributing to evolution of the f/*^'^-' Yukawa matrices are those depicted in 
Fig. iMl The three diagrams with crosses there correspond to contributions of the following 








E<lql<E+dE 

FIG. 14. Diagrammatic representation of the RGE for the up-type Yukawa matrix U. Only the IPI 
scalar exchange [(e)] and the effects of the renormalizations of the Yukawa matrix, of the scalar fields and 
the quark fields in the Yukawa couplings [(b), (c), (d), respectively] contribute when the cutoff is changed 
from E (RHS) to E+dE (LHS). Note that dU stands for U{E+dE)-U{E), etc. The contributions of the 
gauge boson exchanges were not considered in the Figure. 

changes in the Yukawa coupling terms: 

• Yukawa matrix change (renormalization) dlj'^^^ [= U^''\E+dE)-U^^\E)] - Fig. U(b); 

• the scalar field renormalization d(j)^J''> [= (j)^J'\E + dE) — (j)^J'\E)] - Fig. |14|(c); 

• the quark field renormalization du^^^ and du^^^ - Fig. 0(d). 

Figure |1^ is a diagrammatical representation of the physical requirement that the three- 
point (quark-antiquark-scalar) Green function, at one-loop level, in the theory with the 
cutoff E + dE [left side of Fig. |T1|: (a) + . . . -|- (e)] be the same as in the theory with the 
slightly lower cutoff E (right side). 
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Using the results of this and the previous Subappendices, we can then write down the 
one-loop RGE for U^'^^ corresponding to Fig. 0, at the right-handed component [oc(l+75)] 
of the three-point Green function 



(A.23) 



Taking the left-handed component of the three-point Green function results in the Hermitean 
conjugate of Eq. ( |A.23| ), i.e., in an equivalent relation. The first sum on the left (oc A^c) of 
Eq. (pOSi) corresponds to Fig. |T|(c) [cf. Eqs. (|A]Tg)-(^]T|)], the second sum to Fig. ^{d) 
[cf. Eqs. (lAlQl ), (^H], and the third sum to Fig. |l|(e) [cf. Eq. (|A:22|) ]. The left-handed 
part of the Green function yields just the Hermitean conjugate of the above matrix relation. 
The analogous consideration of the three-point Green functions with the down-type external 
quark legs d^^^ and d^^'' gives relations which can be obtained from the above relation again 
by the exchanges f/*-^-* ^ D^^^ and U^'^^'^ ^ D^^^"^ . These relations can be rewritten in a more 
conventional form 



2 e=i 



(A.24) 



and an analogous RGE for D^''\ These RGE's still don't contain one-loop effects of exchanges 
of gauge bosons. However, since the couplings of quarks and the Higgs doublets to the gauge 
bosons are identical to those in the usual MSM, 2HDM(I) and 2HDM(II), their contributions 
on the right of the above RGE's are identical to those in these theories. Therefore, the final 
form of the one-loop RGE's for the Yukawa matrices in the general 2HDM now reads 

dinE [ £=l 

1 2 2 



=1 



(A.25) 



d 



16yr2^^^('^)(E)= AT.VTr ^('^^^Wt + ^W^Wt 
ami? 1 ^ 
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-2^ [t/W^^'^^^W] - AnD^'^A , (A.26) 

where the functions Au and Ad, characterizing the contributions of the gauge boson ex- 
changes, are gauge independent and are the same as in the MSM, 2HDM(I) and 2HDM(II) 

Au-^^-^^^gl + lgl+^^gl, An^Au-gf, (A.27) 
and the gauge couphng parameters gj satisfy the one-loop RGE's 

ISvr^^^. = -C^g] , (A.28) 

with the coefficients Cj being those for the 2HDM's (Nh — ^) 

1 2 1 10 

C3 = -(llTVc - 2n,) , C2 = 7--n, , C,^----n,. (A.29) 

Here, Uq is the number of effective quark flavors - e.g., for E > rrit we have Ug ~ 6; for 
mi,<E <mt we have n^paS, etc. 

This completes the derivation of the one-loop RGE's. 
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